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Abstract 

A stochastic model for the continuous nondemolition ohservation 
of the position of a quantum particle in a potential field and a bo- 
son reservoir is given. lt is shown that any Gaussian wave function 
evolving according to the posterior wave equation with a quadratic 
potential collapses to a Gaussian wave packet given by the stationary 
solution of this equation.. 

The recently developed methods of quantum stochastic calculus [7, ?] can 
serve for the description of the time-development of continuously observed 
quantum systems [1, 4, 5]. 

We apply this approach to describe the time-behaviour of a one-dimen- 
sional quantum particle in the field of the linear force F = kx + mg. The 
effect of coupling of the particle to a measuring apparatus is represented 
by an extra stochastic force. Assuming a non-ideal indirect observation of 
the particle position one can select such an observation channel that the 
observation is nondemolition [4, 5], i.e. it does not affect the actual as well 
as the future states of the perturbed particle. The posterior dynamics of the 
observed particle is then given by the nonlinear stochastic wave equations 
rigorously derived by quantum filtering method in [4, 3]. 

*M.I.E.M., D. Vusovsky 3/12, Moscow 109028, USSR 
ílnstitute of Physics, N. Copernicus University, Torun, Poland 



1 



It is shown that the Gaussian wave packet evolves to the asymptotic sta- 
tionary one the width of which (identified with the Standard deviation) in the 
coordinate representation is given by the formula r q = h/2m[(K 2 + À 2 ) 1 / 2 — 
/t] 1//4 , where k = k/h and À is the accuracy coeficient of the nondemoli- 
tion measurement of the particle position. This phenomenon which cannot 
be explained by the Schròdinger equation (describing the time-evolution of 
the unobserved quantum system) belongs to the class of phenomena called 
watchdog effects [?] . 

Let us assume that the measuring apparatus is modelled by the Bose field. 
The motion of the particle in a potential is distorted by the apparatus so 
the time-derivative P of the momentum of the particle is no longer equal 
to F(X) = —(f)'(X) (0' = d(p/dx). The Heisenberg equations describing the 
time-development of the momentum and the position of the observed particle 
have the form [6] 

P(t)-F(X(t)) = f(t), X(t) = -P(t), (1) 

m 

where the force f(t) in our model is taken as 

f(t) = -{\/2) l/2 [a\t) - a(t)} = (2A) 1/2 ^a t (í), A > 0. (2) 

% 

Here a(t) = b (t), a^(t) = òj(í) are the annihilation and creation quantum 
noise operators with the canonical commutation relations 

[a(t),a\t)} = 5(t-t'), [a(t),a(t')} = 0. (3) 

given by the Standard boson field operators b s , b\, on the half of line s < 
with free evolution b s (t) = b s ^ ct , b\ = b\_ ct at s = 0. 

The position of the particle is assumed to be observed indirectly, together 
with some noise e(t) (error), therefore the measured quantity is 

y(t) = {2X)^X(t)+e(t). (4) 

it is easy to check with the help of (3) that if 

e(t) = a(t) +a f (t) = 23f?a(t) (5) 
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then for P(t) and X(t) satisfying (1) the following commutation relations 
hold, 

[P(t),y(t')}=0, [X(t),y(t')] = 0, Vt'<t. (6) 

From these relations it follows that the preparation for the measurement of 
any functional Y of the past operators y(t'), t' < t does not affect P(t) and 
X(t) as well as any other Heisenberg operator Z(t') of the particle at t' > t. 
In other words, the condition (6) means that the measurement of Y disturbs 
a priori neither the present nor the future state of the particle. 
Note, that y(t) = y\t) and 

b(í),y(í')]=0, Ví,í', (7) 

therefore y(t) can be continuously measured in time as a classical quantity. 
Following refs. [4, 2], one can say that the continual measurement of y is 
nondemolition with respect to the time-evolution of the system. 

Obviously, the minimal distortion of motion (1) will be obtained for the 
Bose field in a vacuum state. In such a situation e and / are white noises, 
e has Standard intensity 1, while the intensity of / is proportional to the 
measurement accuracy: 

</(*)> = <e(0> = 0, (f(t)f(t>))= 1 -\h 2 5(t-t>), 

(e(t)e(t>)) = 5(t-t>). (8) 

It is convenient to rewrite the equations of motion (1) in the form of the Ito 
quantum stochastic differential equations [7, 4, 5, 6]. We obtain 

dP(t) = ^(x(t)),P(t)}dt + (X/2y/ 2 [X(t),P(t)}[dA(t)-dA\t)}, 

dX(t) = ^[P 2 (t)/2m,X(t)}dt, (9) 

where dA(t) — A(t + dt) — A(t) is the stochastic differential of the Standard 
quantum Brownian motion, A(t)\ t=0 = 0, with the generalized derivative 
a(t) = dA(t)/dt. Using the quantum Ito formula [7] one can obtain from 
eqs. (9) the time-evolution equation for any observable Z (a self-adjoint 
polynom in X and P) 

dZ(t) = {-(V^)[^(í),P 2 (í)/2m + 0(x(í))]-ÍA[X(í),[X(í),Z(í)]]}dí 
+ 2ïï{(\/2) 1 / 2 [X(t),Z(t)]dA}. (10) 
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Eq. (10) describes the prior stochastic dynamics of the particle considered as 
an operi quantum system; the term "prior dynamics" means that the process 
Z(t) is not conditioned by the results of the observation. Obviously, the 
prior state of the particle in the Schròdinger picture is described by a mixed 
density matrix even if the initial state is pure (given by the particle wave 
function \1/ and the Fock vacuum vector). 

The Bose field does not only disturb the system but also conveys some 
information about it. This information is contained in the "output field" 
o(í)out — the field after interaction with the system in question. The mea- 
sured quantity y(t) appearing in (4) can be interpreted as 23?a(í) ou t while 
e(t) = 25Ra(í) with a(t) being the input field [6]. Let us rewrite (4) in the 
form of Ito quantum stochastic differential 



for the integral Y(t) = ^y{t')àt'. 

The posterior mean vàlues z(t), i.e. the mean vàlues of the process Z(t) 
which is partially observed by means of Y(t), are defined as conditional 
expectations z(t) = e t (Z(t)) with respect to the observables Y t = {Y(s)\s < 
t. For a given Z, z is a non-anticipating functional of the observed trajectòries 
1 — °f the output process Y. If the initial state of the system is 

pure and the initial state of the bath is a vacuum state then the posterior 
expectation vàlues are realized with the help of the stochastic wave function 
called the posterior wave function. The posterior state is therefore a pure one 
[5]. According to refs. [4, 5], the posterior stochastic wave function <p(t,x) 
satisfies a new nonlinear stochastic (posterior) wave equation which in our 
case has the form 



dY(t) = (2A) 1/2 X(í)dt + 29Í cL4(í) . 



(11) 
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qf \<p àt = (A/2) 



1/2 (x-q)xipdY , (12) 



with the initial condition ip(0,x) 



if)(x). In the above formula 




(13) 



and 



dY(t) = dY(t) - (2\) 1/2 q{t)dt 



(14) 
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denotes the ito differential of the observed commutative innovating pro- 
cess, which is equivalent to the Standard Wiener one as (dY(t)) = and 

(dY(t)f = dt. 

Let us now discuss the problem of the time-development of the posterior 
wave function at t — > oo. We shall assume that the initial state is the 
Gaussian wave packet, 

iP(x) = ( 2 7ra 9 2 ) 1 / 4 exp(- i ^ (x - q) 2 + '-px^j , (15) 

where p and q denote initial mean vàlues of position and momentum of the 
particle and a 2 stands for the initial dispersion of the wave packet in the 
coordinate representation. We shall prové that the solution 0(t,x) of eq. 
(12) corresponding to the initial condition (15) has the form of a Gaussian 
packet, 

<p(t,x) = c(t) exp{(l /h)[rnu(t)[x - q{t)} 2 /2 + ip(t)x]} , (16) 
with the posterior mean vàlues q(t) (given by (13)) and 

h r 

p(t) = - / 0*(t,x)0'{t,x)dx (17) 
i J 

fulfilling linear filtration equations and u(t) satisfying the Riccati differential 
equation. The normalization factor c(t) = (2-kt 2 )^ 1 / 4 up to an inessential 

phase multiplier and r 2 = q 2 — q 2 is a posterior dispersion of position. 

For this purpose it is convenient to rewrite eq. (12) in terms of complex 
osmotic velocity. Let us first introduce 

T(t, x) = R(t, x) + iS(t, x) = h\n <p(t, x) . (18) 

Then by Ito's rule 

dG(0) = G'(0)d0 + ±G"(m0) 2 

applied to the function G = hlmp and taking into account that (dy3) 2 = 
| X(x — q) 2 (p 2 dt we obtain eq. (12) in the from 

dT+[]-h\(x-q) 2 + i<i)-(i/2m)(T' 2 + hT")]dt = {\/2) 1/2 h{x - q)dY . (19) 
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In terms of complex osmotic velocity 

W(t, x) = — T'(t, x) = U(t, x) + iV(t, x) (20) 
m 

eq. (19) can be rewritten as 

dW + [(h\/m){x -q) + (i/h)<p f - i(WW' + hW"/2m)]àt = (X/2) 1 / 2 — dY . 

m 

(21) 

We look for the solution of (21) corresponding to the linear force F(x) = 
Tikx + mg and to the initial condition 

W(0,x) = ^^- = ^(x- q ) + ^ P (22) 
in the linear form 

W(t,x)=w(t)+u(t)x, (23) 
where in accordance with (16) 

w = —ujq H p . (24) 

m 

By inserting W = w, W" = into (21) we obtain the following system of 
equations for the coeficients w(t) = W(t,0) and oo(t) = W'(t,0), 

dw(t) - i[g + u(t)w(t)]dt = (X/2) 1 / 2 — dY(t) , (25) 

m 

with«)(0) = ^q + ip, 

u(t) + — = i[hK/m + iü 2 (t)] , (26) 
dí m 

with uj(0) = — 2 ^ 2 ; which define the solution of (21) in the form given by 
(23). From (24) we get q(t) = -$iw{t) /2$iuj{t) which is the root of the 
equation R'(t,x) = mU(t,x) = for which the maximum of the posterior 
density \0(t, x)\ 2 = exp[(2 / h) R(t, x)] is attained. The posterior momentump 
coincides with mV(t, q) = S'(t, x) x= q and by (24) p(t) = m^s[w(t) +u(t)q(t)}. 
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Eqs. (24)-(26) give the Hamilton-Langevin equations describing the time- 
development of posterior mean vàlues of position and momentum, 

pdt-mdq = h(\/2y/ 2 §;, q(0) = Q, 

(27) 

dp - mgdt = h(nqdt - (A/2) 1 /^ dY) , p(0) = p . 

The posterior position and momentum dispersions for the posterior wave 
function in the form (16) are given by the formulas 

with üu(t) being the solution of eq. (26); therefore the Heisenberg inequality 
r (? 2 (í)r p 2 (í) > h 2 /A is fulfilled. 

The general solution of eq. (26) reads 

w ia + üü(0) tanh(at) ' w ' v ; v ; 

Obviously, lim^oo u(t) = ia, i.e. ia is the asymptotic stationary solution of 
(26). Consequently, posterior dispersions of position and momentum tend to 
finite limits independent of their initial vàlues, from (28) and (29) we get 

Cloc) = — : I — \' 2 \(>< 2 - A 2 ) 1/2 - h] ? 



o(oo) = 7t^=(t^-) K^ + A 2 ) 1 / 2 - 
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1/2 / 2 ,2 X 1/2 

*<«>> = Ü^^J (ï^J ■ w 

Let us pay attention to the particular watchdog effects which can be 
obtained from (30). 

(a) For k < 0, i.e. for the harmònic oscillator we find that the width of 
the stationary Gaussian packet is smaller than that for the unobserved 
oscillator. If the accuracy A — > oo then r 2 (oo) — > but r 2 (oo) — > oo. 

(b) For /t = 0, i.e. for a particle in a homogeneous (gravitation or elèc- 
tric) field, in particular for a free observed particle (g = 0) the Gaus- 
sian packet does not spread out in time. The asymptotic localization 
r 2 (oo) = (h/2m\) 1 / 2 is inversely proportional to the mass of the parti- 
cle and the measurement accuracy coemcient. 
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(c) For k > 0, i.e. for the case of a linear active system (harmònic acceler- 
ator) we obtain a similar watchdog effect as in case (b). 

Thus the collapse problem of the wave function for a quantum particle 
under the position measurement has found the dynamical solution. 
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